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1. INTRODUCTION 

There is widespread agreement that the obUquity and 
precession of the Earth's orbit result in a significant por- 
tion of global temperature variance at frequencies near 
40,000 and 20,000 years. Aside from this "external" forc- 
ing of the climate system, there is no concensus on the 
physical processes governing the remainder of the vari- 
ance. It is generally agreed that the variance result- 
ing from orbital forcing is superimposed on a continuous 
"background" spectrum with a variance that increases 
with decreasing frequency (Kerr 1978). Most attempts to 
model this background spectrum assume a nonlinear re- 
sponse of the climate system to the orbital forcing (Bhat- 
tacharya, Ghil, and Vulis 1982; Imbrie and Imbrie 1980; 
Le Trent and Ghil 1983; Nicolis 1993). Each of these 
studies has succeeded in reproducing only a limited por- 
tion of the spectrum of climatic variations. 

Several studies have recognized that the low-frequency 
portion of the power spectrum, S{f), of temperature vari- 
ations inferred from ice and ocean cores has the form of 
a Lorentzian distribution: (Hasselmann 1971; Komintz 
and Pisias 1979; Lovejoy and Shertzer 1986) 



S{f) 
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This function is constant for very low frequencies and 
proportional to /^^ for / ^ /q. Hasselmann (1971) and 
Komintz and Pisias (1979) have suggested that a stochas- 
tic model may be the most appropriate model for this 
spectrum since a Lorentzian spectrum can result from 
dynamics in which a macroscopic variable is the sum 
of uncorrelated random pulses with a negative feedback 
mechanism limiting the variance at low frequencies. 

At shorter time scales it is recognized that there is 
"persistence" or correlations in time in meteorological 
time series over a range of time scales. Persistence means 
that warm years (or weeks or months) are, more often 
than not, followed by warm years and cold years by cold 
years. Hurst, Black, and Simaika (1965) presented stud- 
ies of these correlations using the rescaled-range tech- 
nique (see Feder 1991 for an introduction). Hurst found 
that time series of annual mean temperature produced a 
power-law rescaled-range plot with an average exponent 
of 0.73. No persistence would yield an exponent of 0.5. 
Numerical studies have shown that series with a Hurst 
exponent of 0.73 have a power spectrum proportional to 
(Higuchi 1990; Gomes da Silva and Turcotte 1994; 
Malamud and Turcotte 1995). 

In this paper we present a model that provides a spe- 
cific physical mechanism for the entire background spec- 



trum from time scales of 1 day to 1 million years. We 
present spectral analyses of paleoclimatic proxy data and 
instrumental data that support the observation of a low- 
frequency Lorentzian spectrum and a higher-frequency 
spectrum proportional to . At very high frequencies 
we found a distinct difference between the power spectra 
of continental and maritime stations. Continental sta- 
tions exhibit power spectra proportional to /~2 at time 
scales less than one month while the spectra of maritime 
stations remain proportional to f~2 down to time scales 
of one day. Our model is based upon an analytic ap- 
proach to modeling the stochastic diffusion of heat in the 
atmosphere and ocean. The difference between conti- 
nental and maritime stations arises because the air mass 
above maritime stations exchanges heat with both the at- 
mosphere above and the ocean below while the air mass 
above continental stations exchanges heat with only the 
atmosphere above it. 

"If turbulent transfer in a system is dominated by ed- 
dies much smaller than the system size, random con- 
vective action of turbulent eddies will be analogous to 
the molecular agitation responsible for molecular diffu- 
sion" (Moffatt 1983). In this approximation, turbulent 
transfer can be modeled as a stochastic diffusion process 
(Csanady 1980). This is equivalent to the Lagrangian 
theory of turbulence for times long compared to the La- 
grangian time scale (time below which particle velocities 
are autocorrelated) (Tennekes and Lumley 1972; Gifford 
1982). In models of climate change, it is common to 
model the turbulent transfer of heat as a deterministic 
diffusion process (Ghil 1983). The stochasticity of tur- 
bulent transfer results in temperature fluctuations from 
equilibrium not present in a deterministic model of turbu- 
lent heat transfer. We will present the power spectrum of 
temperature fluctuations from equlibrium resulting from 
stochastic heat transport in a two-layer geometry appro- 
priate to the atmosphere and ocean. 

The model we present was first solved by van Vliet, 
van der Ziel, and Schmidt (1980) to determine the power 
spectrum of variations due to the stochastic diffusion of 
heat in a metallic film in thermal equilibrium with a sub- 
strate. Temperature variations in the film and substrate 
occur as a result of ffuctuations in the heat transport by 
electrons undergoing Brownian motion. The top of the 
film absorbs and emits blackbody radiation. In this paper 
we use van Vliet et al.'s model exactly as they presented 
it with the atmosphere as the metallic film and the ocean 
as the substrate. Turbulent eddies in the atmosphere and 
ocean are analagous to the electrons undergoing Brown- 
ian motion in a metallic film in contact with a substrate. 
The model studied by van Vliet et al. (1980), with phys- 
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ical constants appropriate to the atmosphere and ocean, 
yield a power spectrum in agreement with that of chmatic 
variations recorded in instrumental records and inferred 
from ocean and ice cores from time scales of one day to 
one million years. 

2. OBSERVATIONS OF CLIMATIC VARIATIONS 
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Figure 1. Logarithm of the power spectrum of atmospheric 
temperature variance predicted by the model as a function 
of the logarithm of the frequency in years" ^. The crossover 
frequencies observed in the climatological record are, from left 

to right, /o — 40,000Hears ' ^ JxtOOyears ' f'^ ^ Imonth ' 

In Figure |l| we present of the logarithm of the power 
spectrum of the variations in temperature as a function 
of the logarithm of the frequency predicted by the model 
and supported by observational data. At low frequencies 
the power spectum is constant. Above / w 40 ooo years 
the power spectrum is proportional to f~^. Above 
/ w 2 oao\ears powcr spcctrum is proportional to 

At very high frequencies (above / w Yn^mFh) 
spectrum varies as for continental stations and re- 
mains proportional to /~2 for maritime stations. The 
crossover frequencies quoted above are those observed in 
the climatological record. The model predicts crossover 
frequencies which are close to those observed in the cli- 
matological record (two out of three are within a factor of 
two of the observed frequencies). All of the power spectra 
that we present in this paper are plotted after taking the 
logarithm of the power spectrum and of the frequency 
against which the spectrum is plotted. All power-law 
functions appear as straight lines with a slope equal to 
the exponent of the power-law. The frequency unit of all 
observational data is years" ^. 
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Figure 2. Logarithm of the normalized Lomb periodogram 
of the temperature inferred from the Vostok ice core as a 
function of the logarithm of the frequency in years"^ 



Figure |^ shows the logarithm of the normalized Lomb 
periodogram of the Vostok SD record extending back 
220,000 years converted into degrees Celsius by the con- 
version factor 5.6% per degree Celsius (Jouzel et al. 
1987). We obtained the spectrum from M. Ghil and 
P. Yiou (1994). It is not possible to directly esti- 
mate the power spectrum of the Vostok record with the 
FFT since the data are unevenly sampled. Numerical 
Recipes suggests the use of the Lomb Periodogram for 
such data (Press et al. 1992). The periodogram shows 
a constant low-frequency region that changes to a re- 
gion proportional to /"^ above / « 40.000 years - Above 
/ w 2 000 \ears power spcctrum of temperature varia- 
tions is oc /~2 . To reduce the scatter of the periodogram 
at high frequencies, we averaged the periodogram in 
logarithmically-spaced bins of size log/ = 0.01 above 
log/ = -4.0. 
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Figure 3. The average unnormalized power spectrum of 94 
monthly temperature time series on a log-log plot as a func- 
tion of frequency in years" ^. 
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In Figure ^ we present the logarithm of the unnor- 
mahzed power spectra of temperature time series as a 
function of the logarithm of the frequency taken instru- 
mentally at higher frequencies. We plot the logarithm 
of the average power spectrum of time series of monthly 
mean temperature from 94 stations worldwide with the 
yearly trend removed. The power spectra were computed 
by taking the modulus squared of the complex Fourier co- 
efficients obtained by computing the Fast Fourier Trans- 
form using the Numerical Recipes routine realft (Press 
et al. 1992). We computed the power spectra of all com- 
plete temperature series of length greater than or equal to 
1024 months from the climatological database compiled 
by Vose et al. (1992). The yearly trend was removed 
by subtracting from each monthly data point the aver- 
age temperature for that month in the 86 year record for 
each station. All of the power spectra were then averaged 
at equal frequency values. The data yield a straight-line 
least-square fit with slope close to -0.5 indicating that 
iS'(/) oc in this frequency region. This is consistent 
with the results of Hurst, Black, and Simaika (1965) who 
analyzed correlations at instrumental time scales with the 
rescaled-range method. Numerical studies have shown 
that a time series which yields a power-law rescaled-range 
plot with Hurst exponent H = 0.73 has a power spectrum 
S{f) oc /~2 (Higuchi 1990; Gomes da Silva and Turcotte 
1994; Malamud and Turcotte 1995). /~2 power spectra 
have also been observed in variations of atmospheric hu- 
midity on time scales of days to years (Vattay and Harnos 
1994) and in tree-ring widths (a proxy for precipitation) 
up to several thousands of years (Pelletier 1995). 
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Figure 4. Averaged unnormalized power spectrum of 50 
continental and 50 maritime daily temperature time series 
on a log- log plot as a function of frequency in years" ^. The 
crossover frequency for the continental specta is /2 = i^p^^fa • 

In Figure ^ we present the average unnormalized power 
spectrum of time series of daily mean temperature (esti- 
mated by taking the average of the maximum and mini- 
mum temperature of each day) from 50 continental and 
50 maritime stations over 4096 days estimated by com- 
puting the Fast Fourier transform as before. Maritime 



stations were sites on small islands far from any large 
land masses. Continental stations, conversely, were well 
inland on large continents, far from any large bodies of 
water. We chose 50 stations at random from the complete 
records (those with greater than 4096 nearly consecutive 
days of data) provided by the Global Daily Summary 
database compiled by the National Climatic Data Cen- 
ter (1994). We found no records in the database that 
were without at least a few weeks worth of missing days 
over the 14 years of data covered by the longest of the sta- 
tion records. We filled in the missing days with the same 
value as the previous day. This will not introduce any dis- 
cernible error since the erroneous datapoints make up at 
most a fraction of one percent of each series. The yearly 
periodicity in these data were removed by subtracting 
from each station's temperature time series a converged 
least-squares fit (using Numerical Recipes (Press et al. 
1992) routine mrqmin) to a function of the form 

T, = + A cos(— i + 0) (2) 
365 

where i is the number of the day in the year and Tav^ 
A, and (j) were the fitting parameters. This procedure 
is a standard one for subtracting the yearly periodicity 
in a meteorological time series (Janosi and Vattay 1992). 

3 

Continental stations exhibit a high-frequency region. 
Maritime stations exhibit /^^ scaling up to the highest 
frequency. 

3. TURBULENT TRANSPORT AS A 
STOCHASTIC DIFFUSION PROCESS 

In most models of climate change, turbulent transfer of 
heat energy in the atmosphere and ocean is modeled by 
a diffusion process (Ghil 1983). This assumption is sup- 
ported by the self-consistent results of experimental stud- 
ies of tracer dispersion which parameterize transport as 
a diffusion process. The resulting measurements of eddy 
diffusivity are nearly constant as a function of depth in 
the ocean and height in the atmosphere. If heat transport 
was dominated by convection currents of the same scale 
as the height of the atmosphere or the depth of the ocean, 
the resulting diffusion coefficient from such large-scale 
tracer studies would be greatly inhomogeneous. Studies 
of the dispersion of Tritium vertically in the ocean have 
resulted in a depth-independent (except for the mixed 
layer extending down to about 100 m) eddy diffusivity of 
about 6x10"^ m^/s (Garrett 1984). Vertical diffusivity 
in the atmosphere is considerably more uncertain, but 
a height-independent order of magnitude estimate of 1 
w? /s for stable air conditions has been quoted by a couple 
of studies (Pleune 1990; Seinfeld 1986). Radar studies of 
trace gases in the middle atmosphere obtained the same 
diffusivity to within an order of magnitude from 5-10 km 
(Fukao et al. 1994). Since the time scale of horizontal 
diffusion in the atmosphere and ocean is so much smaller 
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than the time scale of vertical diffusion, diffusion of heat 
into and out of a local air mass is one-dimensional. For 
this reason, we consider only the variations in local tem- 
perature resulting from heat exchange vertically in the 
atmosphere and ocean. 

Although coherent air motions exist which lead to 
large-scale periodic motions of the atmosphere, these 
motions are strongest horizontally and are predominant 
above the troposphere where less than half of the heat ca- 
pacity of the atmosphere resides (Dunkerton 1993). We 
will assume that such oscillatory air motions transport at 
most a small fraction of the total heat transported ver- 
tically in the atmosphere and ocean. If most of the heat 
is transported by eddies much smaller than the height 
of the atmosphere and the depth of the ocean, turbulent 
transport can be modeled as a diffusion process. 



4. FLUCTUATION THEORY OF CLIMATE 
CHANGE 

A stochastic diffusion process can be studied analyti- 
cally by adding a noise term to the flux of a deterministic 
diffusion equation (van Kampen 1981): 
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J = -o— — h n\x,t) 
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(3) 



(4) 



where AT is the fluctuation in temperature from equi- 
librium and the mean and variance of the noise is given 
by 



(7y(x,t)) =0 



Figure |^ illustrates the geometry of the layer exchanging 
thermal energy with diffusing regions above and below 
it. A diffusion process has a frequency-dependent corre- 
lation length X = (^)h (Voss and Clarke 1976). Two 
different situations arise as a consequence of the length 
scale, 21, of the geometry. For very high frequencies, 
A << 21. In that case, the fluctuations in heat flow across 
the two boundaries are independent. For low frequencies, 
X » 21 and the boundaries fluctuate coherently. First 
we consider high frequencies. Since the boundaries fluc- 
tuate independently, we can consider the flow across one 
boundary only. The flux of heat energy is given by eq. 
(0). Its Fourier transform is given by 



J{k,uj) 



iujrj{k, to) 
ak^ — iuj 



(7) 



where a = — is the vertical thermal diffusivity. The flux 

pc *^ 

of heat energy out of the layer at the boundary a,t x — I 
(the other boundary is located at x = —I) is the rate 
of change of the total energy in the layer E{t): ii^iii — 
J{l,t). The Fourier transform of this equation is 



dt 



E{uj) = - 



{2TT)iuj 



dk e'^^J{k,uj) 



(8) 



Therefore, the power spectrum of variations in E{t), 
Se{oj) =< \E{u;)\'' >, is 



dk 



(9) 



Since AT oc AE, 5't(w) oc lu~^ also. 

If we include the heat flux out of both boundaries, the 
rate of change of energy in the layer will be given by the 



difference in heat flux: 



dE{t) 
dt 



(5) Fourier transform of E{t) is now 



J{l,t) - J{-l,t). The 



{rj{x, t)7j{x',t')) oc a{x){T{x)fS{x - x')S{t - t') (6) 



E(t) 



21 



y I I 

Figure 5. Geometry of the diffusion calculation detailed in 
the text. 

We will calculate the power spectrum of temperature 
fluctuations in a layer of width 21 of an infinite, one- 
dimensional, homogeneous space. The presentation we 
give is similar to that of Voss and Clarke (1976). The 
variations in total heat energy in the layer of width 21 
is determined by the heat flow across the boundaries. 



E{^) = 



(27r)2w J- 



Then, 



S't(w) oc Ssit^) oc 



dk sin(kl)J{k, lu) 



dk sin^{kl) 
D^k^+uj^ 



(10) 



oc w"5(i _ e~''(sin6' + cos 6*)) 



(11) 



where 6 = (^)^ and ujo — D/2P is the frequency where 
the correlation length is equal to the width of the layer. 
When A << 21, the above expression reduces to St{(^) oc 
uj^2. When A >> 21, Srii^) cx tj^a (Voss and Clarke 
1976). 

In Section 2 we presented evidence that continental 
stations exhibit a high-frequency region and mar- 
itime stations exhibit /~2 scaling up to the highest fre- 
quency. This observation can be interpreted in terms of 
the diffusion model presented above. The power spec- 
trum of temperature variations in an air mass exchang- 
ing heat by one-dimensional stochastic diffusion is oc 2 
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if the air mass is bounded by two diffusing regions and 
is oc /~2 if it interacts with only one. The maritime 
stations have a power spectrum up to the highest 
frequency because the air mass above a maritime station 
exchanges heat with both the atmosphere above and the 
ocean below. The fluctuation calculation appropriate for 
maritime stations is one in which the coherent fluctua- 
tions from two boundaries are considered as in the cal- 
culation of the /~2 spectrum. The air mass above con- 
tinental stations exchanges heat energy only with the at- 
mosphere above it. The calculation appropriate for conti- 
nental stations is the one-boundary model which predicts 
the observed spectrum. At low frequencies, hori- 
zontal diffusive heat exchange between continental and 
maritime air masses limits the variance of the continen- 
tal stations. Low-frequency fluctuations of continental 
stations may correspond to fluctuations in the average 
temperature of the atmosphere above the oceans, giving 
a /~2 spectrum. 

At lower frequencies, the atmosphere and ocean 
achieve thermal equilibrium. The variance in temper- 
ature of the atmosphere and ocean is then determined 
by the radiation boundary condition. The fluctuating 
temperature of the atmosphere adds and subtracts heat 
from the atmosphere and ocean through a linear radi- 
ation boundary condition; the heat flux out of the at- 
mosphere is proportional to the temperature of the at- 
mosphere. This results in temperature and irradiance 
variations with a random walk (/~^) spectrum. 

The fluctuating input and output of heat in the re- 
gion will cause large variations from equilibrium. When 
the temperature of the atmosphere and ocean becomes 
larger than the equilibrium temperature, it will radiate, 
on average, more heat than at equilibrium. Conversely, 
when the temperature of the atmosphere and ocean wan- 
ders lower than the equilibrium temperature, less heat is 
radiated. This negative feedback limits the variance at 
low frequencies resulting in a constant power spectrum 
at very low frequencies. 

To show this, we consider a coupled atmosphere-ocean 
model with an atmosphere of uniform density (equal to 
the density at sea level) in contact with an ocean of uni- 
form density. The height of our model atmosphere is the 
scale height of the atmosphere (height at which the pres- 
sure falls by a factor of e from its value at sea level). 
Figure || illustrates the geometry and constants chosen 
(where a is the vertical heat conductivity, p is the den- 
sity, c is the specific heat per unit mass, a is the vertical 
thermal diffusivity, and g is the thermal conductance of 
heat out of the Earth by emission of radiation). Primed 
constants denote values for the ocean. The physical con- 
stants which enter the model are the thermal conduc- 
tance by emission of radiation and the density, specific 
heat, vertical thermal diffusivity, and depth of the ocean 
and atmosphere. The density and specific heat of air 
and water are well-known constants. We chose an ocean 
depth of 4 km and an atmospheric height equal to the 
scale height of 8 km as used by Hoffert, Callegari, and 



Hsieh (1980) in their climate modeling studies. The tur- 
bulent transfer of heat vertically in the ocean by diffu- 
sion and advection is commonly parametrized in studies 
of climate change as a diffusion process with an effec- 
tive thermal diffusivity as we assume in this model (Ghil 
1983). The diffusion coefficients of the atmosphere and 
ocean we employ were discussed in Section 3. 



2 1 2 x=o 

ATMOSPHERE 

a=1.0m7s 
a=1000W/m°K 
p=lkg/m 
c=1000J/kgK 
w, = 8000 km 



OCEAN 

a'=6xlO"''mVs o'=25 W/m°K 
p'=1000kg/m' c'=4200J/kgK 
w, = 4000km 

1 X=W| +w. 

Figure 6. Geometry and constants of the model described 
in the text. 

The equation for temperature fluctuations in space and 
time in the model is eq.(P 

-^(^) — d^^ — " — (^^^ 



dt 



with 



{r,ix,t))=0 



dx 



(13) 



{T]{x,t)r]{x' ,t')) a{x){T{x)fS{x - x')S{t - t') (14) 

The boundary conditions are that there be no heat flow 
out of bottom of the ocean and continuity of temperature 
and heat flux at the atmosphere-ocean boundary: 



'^1 -n 



AT{x 



AT{x ■ 



dAT _ ,dAT 

dx '^=""1" ~ ^ dx ^'^='"t 



(15) 



(16) 



(17) 



At the top of the atmosphere we will impose a black- 
body radiation boundary condition. Most (65%) of 
the energy incident on the Earth is emitted as long- 
wavelength blackbody radiation from the H2O and CO2 
in the atmosphere (Peixoto and Oort 1992). This heat 



5 



emitted from the atmosphere is dependent on the tem- 
perature of the atmosphere at the point of emission ac- 
cording to the Stefan-Boltzmann law. It is common prac- 
tice to assume that temperature variations from equihb- 
rium are small (the global mean temperature has fluctu- 
ated by only about eight degrees Celcius during the last 
glaciation). Within a linear approximation, the emit- 
ted temperature will be proportional to the temperature 
difference from equilibrium (Ghil 1983). The boundary 
condition at the scale height of the atmosphere (which we 
take to be representative of the average elevation where 
radiation is emitted from the atmosphere) is then 



a^U=o=gAr(x-0) 



(18) 



We will use the value g = 1.7 W/m?K as used by Ghil 
(1983) and Harvey and Schneider (1985). 

van Vliet et al. (1980) used Green's functions to 
solve this model. The Green's function of the Laplace- 
transformed diffusion equation is defined by 



iujG{x, x' , iuj) — a{x) 



d^G{x, x' ,iuj) 
dx"^ 



5{x-x') (19) 



where G is governed by the same boundary conditions 
as AT. This equation can be solved by separating G 
into two parts: Ga and Gf, with x < x' and x > a;', re- 
spectively, where Ga and Gb satisfy the homogeneous (no 
forcing) diffusion equation with a jump condition relating 
Ga and Gb'- 



dGg 
dx 



dGb 
dx 



1 



a{x') 



(20) 



The power spectrum of the average temperature in the 
atmosphere in terms of G is given by van Vliet et al. 
(1980) as: 

SAT^Af)oiRe{ / Gi{x,x',iuj)dxdx') (21) 
Jo Jo 



oc Re{ 







Gib{x, x' , iuj)dxdx' 
Gia{x, x' , iLu)dxdx') 



(22) 



where Gi stands for the solution to the differential equa- 
tion for G where the source point is located in the at- 
mosphere. Re denotes the real part of the complex ex- 
pression. Two forms of Gia and Gib are necessary for 
X located above and below x' , respectively, due to the 
discontinuity in the derivative of Gi created by the delta 
function (the jump condition). The solution of Gi which 
satisfies the above differential equation and boundary 
conditions is 

+ cosh — ^ cosh -4 sinh - + — cosh - 23 

L u L Lg L 



and 



where 



Gib = Gia + — sinh( ^ ^ ] 
a L 



(24) 



K = (smh(-) + cosh(-))— smh(- 



+ (cosh( — ) — smh(— )) cosh( — ) (25) 

and L = (^)^ and L' — (^)^. Performing the integra- 
tion van Vliet et al. obtained 

Sat.M) oc i?e(L^(^tanh(^)((^ - 1) 
crL' L cr 

, ,wi. 2qL coshiwil L) — 1 w^^ 

tanh(— ) ^ , h — ) 

^ L ^ a cosh(wi/L) L ' 



gwi ^ wi^ ^ gL 
^ L ^ a 

— tanh(-^) + (l + — tanh(-^)))-i (26) 



( + tanh — tanh — + — 

a L a L L g 



aU —""^ If ' ' '^^ • Lg 
For very low frequencies, 

,,Wi Wi ,,W2. W2 

tanh( — ) « — ,tanh( — ) « — 

cosh(wi/i) — 1 ^ ^ 1 
cosh(u>i/L) ' 2 L2 

Reducing eq. (p6|), 

SATa^ (/) OC CX 



1 + ^ p+f§ 



(27) 



(28) 



(29) 



which is the low-frequency Lorentzian spectrum observed 
in the Vostok data. The crossover frequency as a function 
of the constants chosen for the model is 



fo 



2-K{wiCp + W2C'p'{l + ^)) 2ttwiW2c'p' 



1 



25,000 years 



(30) 



which is within an order of magnitude of the observed 
crossover frequency of the Vostok data, /o = 40 000^ years ■ 



At higher frequencies 
-)« 



tanh( — ) w — ,tanh(— ) « 1 



cosh(wi/i) — 1 \w\ 
cosh{wi/L) 2 L2 



(31) 



(32) 



then 



Sr.M)^li^)H^)H^)^o.f-i (33) 
2 a cp'a' wipcj 
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as observed. The high and low-frequency spectra meet 

at 

/i = _2_(^ )^{— — )S4S( ^ ^ )3 (34) 



10, 000 years ' ' 

which also agrees well with that observed in the Vostok 

data (/i « 2000 years)- 

The time scales of /i and /o correspond to thermal and 
radiative equilibration of the coupled climate system, re- 
spectively. At time scales of 2000 years, the entire at- 
mosphere and ocean are in thermal equilibrium. Besides 
the Vostok data, observational support for this thermal 
equilibration time is provided by comparisons of Antarc- 
tic and Greenland ice cores. Bender ct al. (1994) found 
that temperature variations in Antarctica and Greenland 
are correlated above time scales of 2000 years and uncor- 
rclatcd below it. This suggests that 2000 years represents 
the time scale of global thermal equilibrium. 

Our model makes a specific prediction of the ability 
of the oceans to absorb an increase in atmospheric tem- 
perature resulting from the greenhouse effect. If the rate 
of increase of atmospheric temperature produced by the 
greenhouse effect is less than the vertical heat conductiv- 
ity of the oceans, all of the heat produced by the green- 
house effect can be stored in the oceans for time scales 
up to 2000 years, the thermal equilibration time of the 
atmosphere and ocean. A 1°C per century warming is 
seven orders of magnitude smaller than the vortical heat 
conductivity, indicating that the ocean can easily absorb 
such an increase. A time scale of 2000 years is much 
larger than other values quoted for the delay in atmo- 
spheric warming due to the absorptive capacity of the 
oceans. Most studies have suggested time scales on the 
order of decades (Schneider and Thompson 1981). 

We have applied the same model presented in this pa- 
per to variations in the solar luminosity from time scales 
of minutes to months (Pelletier 1995). A stochastic dif- 
fusion model of the turbulent heat transfer between the 
granulation layer of the sun (modeled as a homogeneous 
thin layer with a radiative boundary condition) and the 
rest of the convection zone (modeled as a homogeneous 
thick layer with thermal and diffusion constants appro- 
priate to the lower convection zone) predicts the same 
spectral form observed in solar irradiance data recorded 
by the ACRIM project and observed in the climate spec- 
tra reported here. The time scales of thermal and radia- 
tive equilibrium of the solar convection zone based upon 
thermal and diffusion constants estimated from mixing- 
length theory match those observed in the ACRIM data. 

5. CONCLUSIONS 

We have presented evidence that the power spectrum 
of atmospheric temperatures exhibits four scaling re- 



gions. We presented a model originally due to van Vlict 
et al. (1980) proposed to study temperature fluctuations 
in a metallic film (atmosphere) supported by a substrate 
(ocean) that matches the observed frequency dependence 
of the power spectrum of temperature fluctuations well. 
The difference between the high-frequency spectra of con- 
tinental and maritime stations may be interpreted in 
terms of the fact that air masses above maritime stations 
interact with the atmosphere above and the ocean below 
while air masses above continental stations interact with 
only the atmosphere above. 
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